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Abstract. We study the surface magnetization of aperiodic Ising quantum chains. 
Using fermion techniques, exact results are obtained in the critical region for 
quasiperiodic sequences generated through an irrational number as well as for the 
automatic binary Thue-Morse sequence and its generalizations modulo p. The surface 
magnetization exponent keeps its Ising value, fis = 1/2, for all the sequences studied. 
The critical amplitude of the surface magnetization depends on the strength of the 
modulation and also on the starting point of the chain along the aperiodic sequence. 



1. Introduction 

Since the discovery of quasicrystals [1] there has been a growing interest in the 
structural and physical properties of quasiperiodic systems (see e.g. [2-5]). Due 
to their situation, which is intermediate between periodic and random ones, these 
systems are particularly interesting from the point of view of their critical behaviour. 

On the Penrose lattice, an approximate renormalization group study of the 
ferromagnetic Ising model [6] suggested a logarithmic singularity in the specific heat, 
i.e. ordinary two-dimensional Ising behaviour. In the same way, random and self- 
avoiding walks, studied via Monte-Carlo techniques, were shown to keep their usual 
asymptotic properties [7] . 

Exact results have been obtained, either for aperiodic quantum spin chains or 
for the corresponding layered classical counterparts ([8-12] and references therein). 
A condition for quasiperiodicity was established [13] in the case of aperiodic chains 
generated through an inflation rule. The one-dimensional lattice is quasiperiodic, i.e. 
it displays a Fourier spectrum with Dirac delta peaks for a set of Q-values which are 
not integral multiples of any single wave-vector, when the characteristic polynomial 
of the substitution has only one root with an absolute value greater than one. 

Mainly bulk critical properties of aperiodic Ising systems have been worked out 
in the last years. Tracy [9] found that, although for quasiperiodically layered Ising 
systems the specific heat displays a logarithmic singularity, this is no longer true for 
some aperiodic sequences which do not satisfy the criterion for quasiperiodicity. The 
low energy finite-size behaviour of the spectrum of quasiperiodic Ising quantum chains 
was found to be in agreement with the prediction of conformal theory by Igloi [10]. The 
bulk critical exponents deduced from the gaps are those of the periodic system. The 
structure of the finite-size Ising spectrum has been carefully studied for the Fibonacci 
chain in [12]. Besides the conformally invariant scaling of the levels at low energy, a 
multifractal scaling is observed for higher ones. 
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In the present work we examine the surface critical behaviour of aperiodic Ising 
quantum chains. The surface magnetization is calculated exactly in the critical region 
making use of fermion techniques [14,15]. The outline of the paper is as follows. In 
section 2 we recall how the surface magnetization of the inhomogeneous Ising quantum 
chain can be obtained. The case of a quasipcriodic chain generated by an irrational 
number is considered in section 3. This is followed by a study of the aperiodic binary 
Thue-Morse sequence in section 4 and of its generalizations modulo p in section 5. 
The results are discussed in section 6. 



2. Surface magnetization of the inhomogeneous Ising quantum chain 

Let us consider a semi-infinite inhomogeneous Ising quantum chain with modulated 
interaction strength A^- and a constant transverse field h=l. The Hamiltonian can be 
written as 

oo 

^ = -2E['-I + VJ^;+i] (2.1) 

where the it's are Pauli spin operators. This quantum chain also corresponds, in the 
extreme anisotropic limit, to a semi- infinite layered two-dimensional classical Ising 
system with modulated interactions on the bonds perpendicular to the surface. 

Using the Jordan- Wigner transformation [16] the Hamiltonian is rewritten as a 
quadratic form in fermion operators cj~, Cj, which may be diagonalized through the 
canonical transformation [14,17] 

a 

where 77^, and 77^ are new fermion operators. and ^'q, are the real normalized 
eigenvectors of two matrices leading to the excitation energies of the diagonal 
Hamiltonian W = X^^, e^{r]^r]^ — 1/2). They are related through 

{^ + B)^^ = eJ>^ (A + B)ij = -Jy- Vi+ij (2.3) 

Below the critical point, in the ordered phase, the ground-state is two-fold degenerate 
and the lowest excitation Cq, which corresponds to an eigenvector localized near to 
the surface, vanishes exponentially with the size of the system. This provides the 
condition for criticality [18] 



N 

hm n(^)c = l (2-4) 

JV—>oo 



The surface magnetization can be obtained using a simple method due to Peschel [15]. 
The asymptotic form of the surface two-point correlation function gives the square of 
the surface magnetization which can be expressed as the matrix element (Ojcrf jcr) 
where |0) is the vacuum of Ti and [cr) =?7q"|0) its first excited state. Rewriting in 
terms of diagonal fermions through (2.2), one obtains = ^'□(1). Using (2.3) with 
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a = in the ordered phase where the lowest excitation vanishes, a recursion equation 
for the successive components of $o is obtained leading to 



^o(i) = (-ir'^o(i)nV' i>i 



(2.5) 



i=l 



The surface magnetization then follows from the normalization of the eigenvector and 
reads 



oc ] 



-1/2 



(2.6) 



3. Quasiperiodic modulation generated by an irrational number 

Let the position Xj of the j-th spin along the chain be defined as [10, 13] 

fj + l" 



Xi = jb + 



{a -b) (a > 6) 



(3.1) 



where [x] is the integer part of x. When uj is irrational, one obtains a succession of 
short and long bonds which are distributed according to a quasiperiodic sequence. 
With a modulation amplitude rf such that o-bonds have a strength A and 6-bonds a 
strength Ar, the number of modified bonds from the surface to site j is given by 



n-i =J 



J + 1 



UJ 



(3.2) 



corresponding to an asymptotic density p^ = l — l/w. The critical coupling A^ follows 
from (2.4) with 



N 



lim n(^)c= lim A^^''~=l 



SO that 



JV— »oo 



The surface magnetization (2.6) involves the sum of products 



CO J oo oo ^ , s 2j 



j=l i=l 



(3.3) 



(3.4) 



(3.5) 



t In the following we tcike r > without loss of generality. 
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Figure 1 . Surface magnetization rus of the Fibonacci Ising quantum chain with a 
modulation amplitude r = 1, 2, 3, 5 from top to bottom. 



In the last term n -, given by (3.2), was rewritten as 



(3.6) 



where {x} is the fractional part of x. The eigenvector $q, which is localized near to 
the surface in the ordered phase, extends into the bulk in the critical regime (A— + A^). 
Then the last factor in (3.5) may be replaced by its averaged value [11]. Since w is 
irrational, {{j + Vj/w} is uniformly distributed over [0, 1] and the average is given by 



'2x 



- 1 



^2 Yt/x j-2 



(3.7) 



Summing the geometrical series, one obtains the exact leading contribution in the 
scaling region 



,(1) 



1 

a: 



In 5 



2 ll/2 



1 



1/2 



(3.8) 



The surface exponent (3^ = 1/2 is the same as for an homogeneous Ising system and the 
critical amplitude keeps the same form for any sequence generated by an irrational. 
It only depends on the modulation amplitude and the critical coupling. 

When iv is the golden mean (•\/5 + l)/2, the modulation is distributed according 
to the Fibonacci sequence. Such a sequence may be obtained in various ways. It can 
be represented by an infinite word = 01001010... where stands for an a-bond 
and 1 for a 6-bond. This word is the limit of an infinite sequence of finite words 
obtained through the recursion relation 



n+l 



W W 



(3.9) 
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Figure 2. Ratio of the Fibonacci surface magnetization rus to its leading critical 
term m^^^ ~ e^/^ for r = 2, 3, 4, 5 from bottom to top. The linear behaviour 
corresponds to a correction to scaling of order e^/^. 



to which is associated the characteristic polynomial — x — 1 [9]. The length of 
succesive words is given by Fibonacci numbers. 

The Fibonacci sequence also follows from the inflation rule [9] 



^ 01 1^0 (3.10) 

giving 





01 
010 

(3.11) 

01001 
01001010 



The surface magnetization of the Fibonacci quasilattice as given by (2.6) is plotted 

in figure 1 for different vahics of the modulation amplitude r. The ratio m^/rrSp shown 
in figure 2 as a function of e = 1 — A^/ A indicates that the corrections to scaling are 
of order e^/^ in this case. 



4. The binary Thue-Morse sequence 

Let us now consider the binary Thue-Morse sequence /j = or 1 (j = 1,2,3...) which 
is not quasiperiodic [19] and sometimes called automatic since it can be generated 
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Figure 3. Surface magnetization of the binary Thue-Morse aperiodic chain with a 
modulation amplitude r = 1, 2, 3, 5 from top to bottom. 



through a 2-automaton. It may be deduced from the binary representation of non- 
negative integers [20] 

1 10 11 100 101 110 111 ... (4.1) 
by counting the sum of their digits modulo 2 

01101001... (4.2) 
or through the recursion 

W„+i = W„ {W-\ W^) W„_, (4.3) 
It also results from the substitution 

^ 01 1 ^ 10 (4.4) 

leading to 


01 

0110 

(4.5) 

01101001 

0110100110010110 




The aperiodic chain is constructed through the correspondences — > A and 1 — > Ar. 
It is clear from (4.4) that the asymptotic density of modified bonds is = 1/2 and, 
according to (3.3), 

A, = r-V2 (4.6) 
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Figure 4. Ratio of the binary Thue-Morse surface magnetization rus to its leading 
critical contribution m's^ for a modulation amplitude r = 2, 3, 4, 5 from top to 
bottom. The corrections to scaling, due to the higher regularity of the sequence, 
are much weaker than for the Fibonacci sequence. 



This value of the density is also reached for any finite sequence with an even length 
so that the number of modified bonds satisfies 

«2fe = (4-7) 
The substitution (4.4) is such that the binary Thue-Morse sequence is reproduced by 
selecting the sequence of its odd elements {f2k^i = fk+i) indicated in the last line 
of (4.5). As a consequence 

%fc+i=fc + /i+\ fc = 0,1,2... (4.8) 
The sum involved in the surface magnetization (2.6) can then be written as follows 

ooj ooy,s4fc °°/\ \ 2(2fc+l) 

j=li=l k=l ^ ^ fe=0 ^ ^ 

Near the critical point the last term may be replaced, as above, by its average over 
the infinite sequence in which and 1 occur with the same probability. Then 



1-2/, 



(2) \ 1 



2 

and the exact leading critical behaviour reads 

1/2 



= -(r + r-i) (4.10) 



.(1) - 



r ■ 1 



1 - 



(4.11) 



The surface magnetization obtained numerically using (2.6) is plotted in figure 3 for 

different values of r and the ratio m^lrrS}'^ is shown in figure 4. As a result of the 
higher regularity of the Thue-Morse sequence, the corrections to scaling are much 
weaker than for the Fibonacci sequence. 



8 



L Turban and B Berche 



5. Thue-Morse sequences modulo p 

The results of the last section can be generalized by considering Thue-Morse sequences 
fj^^ defined modulo p. These sequences are generated by taking the sum modulo p 
of the digits of non-negative integers written in base p. With p = 3, for example, one 
obtains 

1 2 10 11 12 20 21 22 100 101 102 ... (5.1) 

and the sequence fj'''^ reads 

012120201120... (5.2) 
fj^^ also follows from the substitution 

0^012. ..p-1 1^12. ..p-10 ... p-l~^p-10...p-2 (5.3) 
which for p = 3 leads to 



012 

012120201 (5.4) 
012120201120201012201012120 

One associates a modified coupling on the chain with any nonzero digit along the 
sequence, i.e. we use the correspondence 

O^A 1,2, ...,p-l^ Ar (5.5) 

Then, according to (5.3), the asymptotic density of modified bonds is p^ = {p—\)/p 
and the critical coupling Ag = r~^^^^/^. The asymptotic density is also obtained when 
the number of bonds is any multiple pk of p so that 

%k = pkpoo = Hp - 1) (5-6) 

One may check on the last line of (5.4) that /p^+i = fl^i whereas fpk'+2^fpk+3--- 
are deduced from f£^i through circular permutations on 0, 1, 2, . . . ,p— 1. Then, the 
number of modified bonds after pk + q steps can be written as 

npk+q = Hp - 1) + ™fe+i il) < 5 < p - 1 (5.7) 



where 

9-1 



. — T |_ ^ mod p 



0<q<p-l, mfe(0) = (5.8) 
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Figure 5. Surface magnetization of binary Thue-Morse aperiodic chains with r = 2, 

starting on the j-th digit along the sequence witli j = 1, 2, 3, 4. Tlie amplitude of the 
leading term given by (4.11) when i is odd, is the geometric mean of the amplitudes 
obtained for even values of i. In the strong coupling regime the local environment 
governs the behaviour of the surface magnetization. 



and 

9{i) = 1 - i5i,o 

The sum in (2.6) may be splitted into p parts and rewritten as 



(5.9) 



^ ^ 2pk p-1 ^ 2g 



oo J o 

EnA.-'-E(^j +E(f 

j=li=l fc=l ^ / g=l ^ 



.2</(p- 



i)/p V ( ^ 



2pk 



When A^A^, r-2™'=+i(9) in the last sum may be replaced by its averaged value. On 
the infinite sequence, m^(g) is equal to q—1 with probability q/p and equal to q with 
probability {p—q)/p, so that 



-2mfe+i(g)\ ^ 9^ 



-2(g-l) 



+ ip- q)r' 



-2q 



I 

In the critical region, one finally obtains 

' j,i/p _ j.-1/p 



(5.11) 



,(1) = 



VP 



r — r 



2p- 



1/2 



(5.12) 



for the leading contribution to the surface magnetization in agreement with (4.11) 
when p = 2. 
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Figure 6. Ratio of the surface magnetization tUs to its leading term m^^^ for the 
binary Thue-morse Ising chain with r = 2 and different starting points i along the 
sequence. Although the leading amplitude is the same for all chains with odd values 
of i, the corrections to scaling are different. 



6. Discussion 

In the cases studied so far the modulation was distributed on the spin chain according 
to a given sequence, starting with the first digit of the sequence. One expects a 
change in the surface magnetization amphtude for chains starting elsewhere along the 
sequence. The previous calculations are easily extended to the case where one begins 
the chains on the i-th digit. 

For the Fibonacci sequence, n • in (3.6) is changed into 



leading to 



J_±l 

CO 



Inr^ 


1/2 




_r2 _ I 







1/2 



(6.1) 



i > 1 



(6.2) 



(2) 

With the binary Thue- Morse sequence starting on i even, i.e. with f^i^ one obtains 

n^k = k-l + fi%+fi? n2k+i=k + fS^ 1 = 21 (6.3) 
and then, on the surface site, 

211/2 

r+1 ^ \ X 



mW(2Z) 



(6.4) 
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Otherwise, when i is odd, 

= ^ n^k+i = + fkli i = 2l + l (6.5) 

like in (4.7-8). Then the leading term is still given by (4.11) which is just the 
geometrical mean of the two allowed values in (6.4). The dependence of the surface 
magnetization on i is shown in figure 5. The local environment clearly governs the 
behaviour when X^/X 0. The surface magnetization is reduced when the surface 
site is linked to the bulk via a weak bond for i = 1 and 4. Near the critical point, on 
the contrary, only the long-range structure of the modulation and the parity of the 
surface site arc; important. 

Although the amplitude of the leading contribution to the surface magnetization is 
the same for all the aperiodic chains starting with odd values of i along the Thue- Morse 
sequence the corrections to scaling differ as shown on figure 6. 

Similar results would be obtained with the Thue-Morse sequence defined modulo p 
with the same leading amplitude for chains starting oni = pl + 1 along the sequence. 

For all the sequences considered here, the aperiodic modulation, either quasiperi- 
odic or aTitomatic, leaves the surface magnetic exponent unaffected, /3g = 1/2, like in 
the homogeneous semi-infinite Ising model. On the other hand, the critical amplitude 
depends on the modulation and it would be of interest to study other surface quanti- 
ties in order to check whether surface critical amplitude ratios remain universal. The 
study of aperiodic sequences which do not lead to ordinary bulk Ising behaviour [9] 
would also be worthwhile. 
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